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Abstract 

In this paper we focus on analytical calculations involving null geodesies in some spherically 
symmetric spacetimes. We use Weierstrass elliptic functions to fully describe null geodesies 
in Schwarzschild spacetime and to derive analytical formulae connecting the values of radial 
distance at different points along the geodesic. We then study the properties of light triangles 
in Schwarzschild spacetime and give the expansion of the deflection angle to the second order 
in both M/ro and M/h where M is the mass of the black hole, ro the distance of closest 
approach of the light ray and h the impact parameter. We also use the Weierstrass function 
formalism to analyze other more exotic cases such as Reissner-Nordstr0m null geodesies and 
Schwarzschild null geodesies in 4 and 6 spatial dimensions. Finally we apply Weierstrass 
functions to describe the null geodesies in the Ellis wormhole spacetime and give an analytic 
expansion of the deflection angle in M/h. 



1 Introduction 



Geodesies in Sehwarzsehild spaeetime have been studied for a long time and the importance of a 
good understanding of their behavior is clear. In this paper we shall focus on analytical calcula- 
tions involving null geodesies. While these are interesting in their own right, calculations like this 
are also important for experiments testing General Relativity to high levels of accuracy. Examples 
of two such proposed experiments are " The Laser Astrometric Test of Relativity" or LATOR and 
"Beyond Einstein Advanced Coherent Optical Network" or BEACON, which are both using paths 
of light rays to verify General Relativity and are described in detail in [T] . Both are intended to 
measure second order effects in light bending Elliptic functions have been used to describe the 
geodesies in Sehwarzsehild spaeetime before, mainly in W and more recently in [31[TS]. In pi [T3] 
the focus is mainly on the paths of massive particles and even though they mention the possibility 
of using Weierstrass functions in the null case, they don't go into much detail. The discussion in 
[3] is concerned with null geodesies around a charged neutron star using the Reissner-Nordstr0m 
metric, a case we also study but with a different emphasis. In this paper we begin by providing 
a complete description of Sehwarzsehild null geodesies in terms of Weierstrass functions and then, 
and this is our principal innovation, using various "addition formulae" for Weierstrass functions 
[S], we derive some analytical formulae connecting values of radial distance at different points 
along the geodesic. The motivation is to develop, as far as is possible, optical trigonometry in 
the presence of a gravitating object such as a star or a black hole. To that end we use these the 
additon formulae to study the properties of light triangles in the Sehwarzsehild metric and obtain 
the deflection angle of the scattering geodesies to second order in both M/tq and M/b where M 
is the mass of the black hole, tq the distance of closest approach of the light ray and b the impact 
parameter. 

In the final section we show how the same methods to treat null geodesies in more exotic space- 
times; charged black hole, the Ellis wormhole '6^ and Sehwarzsehild black holes in 4 and 6 spatial 
dimensions. Although not a primary concern of the present paper, it is worth remarking that the 
addition formulae for Weierstrass functions that we make use of are closely related to the existence 
of an abelian group multiplication law on any elliptic curve [7] and suggest, in view of the impor- 
tance of the complex black hole spacetimes at the quantum level, that it might prove fruitful to 
explore this aspect of the theory further. 

The organization of the paper is as follows. In section 2 of this paper, we provide the full 
solution for Sehwarzsehild null geodesies in terms of Weierstrass elliptic functions and apply it to 
obtain addition formulae connecting three points on the geodesic. We then calculate the deflection 
angle of the scattering geodesies to second order in both M/tq and M/b where M is the mass of 
the black hole, tq the distance of closest approach of the light ray and 6 the impact parameter. The 
section is concluded with the discussion of the light triangles and Gauss-Bonnet theorem. In section 
3, we apply the Weierstrass function formalism to further examples such as Reissner-Nordstr0m 
null geodesies and Sehwarzsehild geodesies in more spatial dimensions. At the end of the section, 
we give a detailed description of the Ellis wormhole null geodesies. 

2 Sehwarzsehild null geodesies 

The equation obeyed by a null geodesic r{(j)) in the Sehwarzsehild metric is 



where P — E'^/L'^ — 1/6^. Here E is the energy of the light, L the angular momentum and 
b the impact parameter. Interestingly the same equation arises for a null geodesic r(0) in the 
Schwarzschild-de-Sitter or Kottler metric [5] [5] and many of our results remain valid in that case. 




(1) 
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Geometrically, one may regard solutions of ([T]) as unparameterised geodesies of the optical metric 



dsl = + t4m (de' + Od^^) , (2) 

with 9 — ^. Introducting the isotropic coordinate p — i(r — Af) + ^y/r{r — 2Af), we find that 

dsl = n^{p)\^dp'^ + p^{de^ + sin^ edc/)^)^ , (3) 



where 



(1 ' 
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^(P) - V, ' 't\ ■ (4) 



Thus our results also apply to light rays moving in an isotropic but inhomogeneous optical medium 
in flat space with refractive index n{p). 

Another interpretation of ([T]), recently exploited in ITUj, is provided by substituting r = i in 
and differentiating to obtain 

with 

F{u) = iMh^u^ . (6) 

Now ([5|) is the equation governing the motion of a non-relativistic particle of angular momentum 
per unit mass h moving under the influence of a central force F(u). In our case the effective 
force F{u) is attractive, and varies inversely as the fourth power of the distance. A search of the 
voluminous nineteenth century literature on such problems reveals that it was comparatively well 
known that although this problem admits some simple exact solutions, which we shall detail below, 
the general solution requires elliptic functions. 

If we had adopted isotropic coordinates and substituted p = we would have obtained a very 
different formula for F{u). In fact in that case we would have 



,(l+^)5(l_ M«) 
(1 2 



F{u)^2Mu' ' (7) 



As we shall see in detail in a later section, the null geodsics of neutral Tangherlini black holes 
in D spacetime dimensions correspond, in Schwarzschild coordinates, to the motion of a non- 
relativistic particle with a force F{u) cx . The cases D = 4,5,7 are the only cases known to 
be integrable in terms of elliptic functions. In fact the cases D = A and D = 7 may be related 
by a conformal mapping introduced in his context by Bohlin |12j and elaborated upon by Arnold 
[13] . The Bohlin- Arnold mapping is a type of duality, i.e. it is involutive, and the case Z? = 5 is 
self-dual. 



2.1 Weierstrass functions solution 

Substituting y = M/2r - 1/12 in (P) gives 



where 



{yr = ^y'-^y-93, (8) 



1 /M^ ^ 



^^=216- T)^- 



In the case 173 ^ ±1/216 the general solution to this equation is y{(p) = p{(f) + C) where p{z) is 
Weierstrass elliptic function and C — const. Detailed description of these functions together with 
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the proof of the above statement can be found in [S] pages 429-444 and page 484. For the critical 
values of the equation for r can be integrated to give 



r{(f>) = Mil + cos(0)) (10) 

in the case (73 = 1/216 and 

Mil 

(11) 



r{(j)) 3 1 ± cosh((/)) ' 

in the case (73 = —1/216. The former, geometrically a cardioid in (r, (/>) coordinates, starts at 
the singularity, reaches the horizon from below and then returns back. The latter describes two 
types of trajectories, one starting at infinity, the other at the singularity and both approaching the 
photon sphere, never reaching it. 

Now suppose that 173 ^ ±1/216 and M^P < 1/27. Then the polynomial Ay'^ — y/12 — has 3 real 
roots ei > 62 > 63 and the half-periods of the corresponding Weierstrass function p are 

r dt 

/ =, 12 

4 v/4t3-i/12-g3' 

W3 = -j / = , 13 

i-00 + t/12 - 4*3 

where oji G R and iu)-^ € K. In this case p[z) is real on a rectangular grid with vertices 
0, wi, CJ3, oji + aj3 and since is real, the only physical solutions to (2) are = p(0 -I- 0o) or 
y(0) = p((/) + 00 + W3) where (/)o G M. We have two cases: 

(i) Scattering paths 

The point r = 00 corresponds to y = —1/12 and p(z) takes value —1/12 at z such that Im(z) = a;3. 
Therefore, choosing line (/> to be the axis of symmetry, we have y{(\)) — p(0 + oji + (^3) and so 

M 1 

^-- + 2p(<^-fwi-|-c^3), (14) 
where of course the function p depends on P . Here the range of is [— /3, /?] where 

ie,, V4t^-i/12-g3 
In this notation the angle of deflection 6(f) is Scj) — 2/3 — ir. 

(ii) Trapped paths 

These begin and end at the singularity and r = corresponds to ?; = 00. So, choosing line (j) — 
to be the axis of symmetry once again gives y{4>) — p{(t> + lji) where (p G [— cji, wi]. So in this case 

M 1 

- + 2p((/)-t-wi). (16) 



r((/)) 6 



Figure 1 shows the argument of p in the scattering and trapped cases. 
Now suppose that M^P > 1/27. Then wc have 

(iii) Absorbed paths 

These go from infinity to r = or from r to infinity and therefore the solution is uniquely de- 
termined by P. There is only one real root of the r.h.s. of Weierstrass equation, ei < —1/12 and uji 
defined as before is again a half-period. For each P there is a solution of the form y(0) — p((/) + 0o). 



3 



r = 00 



"3 




Re axis 



r = 



Figure 1 The argument of p in the complex plane, S corresponds to scattering trajectories, T to trapped ones 
^0 € and so by uniqueness, all physical solutions arc of this form. 

We can take = which means defining the line = by the direction in which the path 
leaves/hits r = 0. Then the range of (p is [—a, a] where 



J-i 



dt 



and the solution is 



/i2 V4t3 - t/12 - ga ' 



M 1 „ 

T + 2p(0) . 



r{(j)) 6 



(17) 



(18) 



A diagram of the complex plane corresponding to absorbed trajectories may be found in the Figure 
2. 
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Figure 2 The argument of p in the complex plane corresponding to absorbed trajectories 
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2.2 Addition formulae 



As shown in [5] page 440, Weierstrass functions satisfy an addition formula of the form 



p{x + V) = ^ 



p'{x) - p'{y) 
p{x) - p{y) 



p{x)- p{y) = F{p{x),p{y)), 



where 



F{x,y)^- 



^4x3 - x/l2 - 33 - V4y3 - y/12 - 53 



x-y 



x-y. 



(19) 



(20) 



We can apply this result to null geodesies to obtain an expression for r(0i +02) as a function of 
r{(j)i) and r{4'2)- Of course, if any such formula is to be useful in some experimental setup, we need 
to be able to easily find the line = 0. Also, because of the additive constant in the argument 
of the Weierstrass function, we cannot apply the addition formula directly because the sum of the 
two arguments will not correspond to the sum of the two angles. Fortunately, in the case of the 
scattering and trapped orbits, choosing the line = to be the axis of symmetry takes care of both 
problems. Take the scattering orbit for example. The axis of symmetry is easy to find, and we can 
apply the addition formula for p to 3 points on the orbit j/i = 2/(</>i), 2/2 = y{4>2) and 2/3 = y{(f>3) as 



yC^4'^) = p(^<Pt + 3wi + 3a;3) = i^(F(2/i, 2/2), 2/3) 



(21) 



which works because 2tJi and 2tj3 are periods of p. 

Now, letting ^3 = gives 2/3 = £2 with 62 directly related to the distance of closest approach dmin 
as 62 = M/2dinin — 1/12. Then we obtain an addition formula for 3 points on the orbit in the form 



M 



1 



2r{(j)i + cj)2) 12 



F F 



( M 



1 



M 



V2r(0i) 12'2r(02) 



1 

12 



,62 



(22) 



The same procedure for trapped orbits gives the same formula only with ei instead of 62 where ei 
is related to the maximal attained distance dmax by ei = Af/2(iinax ~ 1/12. 

In the absorbed case, the lack of additive constant in the argument of the Weierstrass function 
means that we can apply the addition formula directly to obtain algebraically simpler result 



M 



2r(0i+02) 12 



M 



1 



M 



2r( 



12 2r( 



1 

12 



(23) 



In this case we can use the euclidean angle between the direction of the ray and 0-direction ^ 
which satisfies 

^ dr ,„ , , 

tan?/'=-— (24) 



Then the addition formula can be written as 
2 



1 



1 



?-(0l) ?'('/'2) 



2M( 



1 



1 



1 



r(0i) r{(j)2) r{(j)i+(l)2) 



tan'!/)(0i) tan'!/)(02) 



r(02) 



(25) 

However, in this case this is not very useful since it is practically impossible to identify the line 
= for such a choice. We could of course make a different choice, like r(0 = 0) = R for some 
chosen R, but then the obtained addition formula would not be analytic anymore because we would 
need to find the corresponding additive constant (j)Q given by the integral 



»o = 



dt 



(26) 
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2.3 The deflection angle 

We start from the equation for u ~ 2M/r which is 



^)^.= -„^.4M=P..= -„^-.3 + .^ ,27, 

del) J 

where /i = 2M /tq, tq is the distance of closest approach, for scattering orbits. Then the deflection 
angle Scf) is given by 

50 = 2/ =-7T = 2I-n. (28) 

The integral can be rewritten in terms of a; — u/ fi which gives 

/(m) = /' '^'^ . (29) 

Jo ^/(^^^^y^Tl^^ 

This integral can be expanded in the powers of fi, for /i sufficiently small, as 

We are only interested in small values of /i and this expansion clearly converges at least for < 2/3 
since (1 — x^)/{l — x^) < 3/2 for x G (0, 1). Calculating the first three terms in this expansion 
results in an expansion for the defiection angle is (substituting for /i) 

4Af /57r 4\ M2 „^ 
50=- + 3(---j^ + O(.-). (31) 

Now, we have P — 1/6^ and so define — 2M/b. Then 

iy2 ^ ^2 _ ^3 ^ (32) 

and working to the second order gives 

fi = iy+^fi'^ . (33) 
Substituting into the expansion for 54> then gives 

<5</. = 2^ + — + 0(^3) = _ + + 0(^,-3) _ (34) 

which is the expansion of the defiection angle to second order in 1/6. 
2.4 Angular sum in light triangles 

Because the Gauss curvature of the optical metric restricted to the equatorial plane is negative, 
the angular sum of a triangle made up of geodesies must less than tt unless the triangle encloses 
the horizon [T3]. One might hope to get a more precise statement using the addition formulae. To 
this end, let be the physical angle between the direction of the light and the (j> - direction. Then 



1 1 dr Pr^-r + 2M , , 

tan9 = —= = \ -VT ■ (35) 

iF rd(^\r~2M ^ ' 



Note that this formula is valid for the Schwarzschild solution but not the Kottler solution with 
non- vanishing cosmological constant [S]. 
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Figure 3 Light triangles 

Now consider 3 light rays, forming a triangle around the origin with Pi, P2 and P3 and vertices 
at the radial coordinate Ri, R2 and R3. To simplify the notation, define 



PiR^^-Rj+2M ^ ^ 

Wij = \ ' ' ■ (36) 



Rj - 2M 



Then from the Figure 3 it is clear that 



a = 7r-tan"^ 14^11 -tan"^ 14^21, (37) 
/3 = TT - tan"^ 1^22 - tan-i W2.2 , (38) 
7 = TT - tan"^ Wi3 - tan"^ W33 • (39) 

Alternatively, if the origin is not inside of the triangle, then from the Figure 3 it follows that 

a = tan-i Wn + tan'^ W21 , (40) 
/3 = TT- tan-i VK22 - tan-i W32, (41) 
7 = tan"^VFi3+tan"^W^33. (42) 

Further analytical work in this general case doesn't seem to lead anywhere, because the distances 
i?2 and i?3 are not independent, but finding a formula for the relation between them is impos- 
sible. We can however consider a symmetric case with all R's and P's equal. Then its angles are 
given by 

1 jPR^-R + 2M 

a = 7r-2tan-iW — . (43 

\ R-2M ^ ' 
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2.5 Gauss-Bonnet theorem 



One of the way to calculate this is transform the optical metric into the form 

ds^ = V + C{pfd4? . 



KdA = 



KdA 

a/2 
-a/2 



An alternative approach to finding the angular deflection is using the Gauss-Bonnet theorem |14) . 
Consider the setup in the Figure 4. Then by the Gauss-Bonnet theorem we have 

a + TT + / KdA = 2n . (44) 

(45) 
(46) 

(47) 
(48) 
(49) 

which holds for any scattering path. It doesn't seem to be very useful when it comes to evaluating 
a but it is an interesting expression. Rewriting this in terms of r gives another interesting identity 

r — M TT 

- rlr — — . 

2Mr 2 



Then 



and so 



Now, 



dC 
dp 



dC 
dp 



r—r{tp) 



dC _dC dr _ r-M 
dp dr dp - 2M 



Therefore we get 



«/2 



r{(f>) ~ M 



-.d4) 



^yjr - 2My/Pr 
where tq is the distance of closest approach. 



(50) 



3 Further applications of Weierstrass functions 
3.1 Reissner Nordstr0m null geodesies 

As mentioned in the introduction, these have been studied previously using Weierstrass functions 
in [3] . In this case, the relevant equation for u = 1/r is 

=P-u^ + 2Mu''-Q^u\ (51) 

One may verify using the formulae in [5] or directly, that just as in the case of the Schwarzschild-de- 
Sitter metrics, so with the Reissner Nordstr0m metrics that the cosmological constant does occur 
in ([FT]) . Thus some of the results in [TT] , which appeared on the archive subsequently to the first 
version of this paper, follow directly from the work of the present section. 

The r.h.s. of (|5ip always has a real root so let xo to be one. Then define s — u~ xq. This gives 



where 



{s'f ^As + Bs^ + Cs^ + Ds^ , 

A = QxlM - 2xq - AxlQ^ , 
B = 6xoM - 1 - 6xlQ^ , 
C = 2M - AxoQ^ , 
D = -Q\ 



(52) 



(53) 
(54) 
(55) 
(56) 
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Figure 4 Scattering light ray 



Now substitution i/j = 1/s takes it into the form 

{ilj'f = Aij;^ + + Cij; + D, (57) 
and finally setting ^j} = Ay/A — B/'iA gives 

{y'f = 4y3 _ g^y - 53 , (58) 

where 

.2 = ^-^, (59) 
ABC A^D , , 

Therefore this time the solution will be given by 

1 ZA 
;M='°+12p(^ + &)-i?' ^''^ 

where is a complex constant. However, the more complicated relation between r and p and also 
many different constants make it algebraically very challenging to analyze the situation any further 
and find a suitable or addition formula similar to the Schwarzschild case. For that purpose, 
consider the equation for r(A) where A is an affine parameter of the path. This equation is 

Now the motion is only possible in regions where /(r) > 0. If these split into two disconnected 
ones, then that must be the case in which the r.h.s. of the corresponding Weierstrass equation has 
3 real roots and we have scattering and trapped paths. If there is only one such region, we know 
we have the case where the above mentioned r.h.s. has only 1 real root and we have absorbed 
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paths. 

First, the roots of f'{r) are 



3M± 

r± = ^ . (63) 



Physically we want > and so r± G M+. Clearly the regions where /(r) > will be 

disconnected (and there will be 2) if P > and f{r+) < 0. Suppose that this is the case and 
consider scattering paths, let a be the distance of closest approach. Then a = Cj and we can write 
y(^) = p{(l) + uji) for some i G {1,2,3} because the path is symmetric. As before, this corresponds 
to choosing the line (f) = to he the axis of symmetry. We can compute a as the largest root of 
/(r) =0 and so we obtain an addition formula 




r{(l>i+(l>2) 



For trapped orbits the same addition formula applies, only in that case a is the largest attained 
distance and is given by the second largest root of /(r) = 0. 

If /(r+) > then we have orbits that go in, miss the singularity, and continue to another asymp- 
totically flat region of spacetime. These satisfy the same addition formula like the scattering ones. 
In the case P = the equation can be integrated and the solution is 



= arctan((^ — <Po) , (64) 



^2Mr -Q'^-r'^ 

where 0o is the constant of integration. Another special solutions solutions can be found when 



3r1 



(65) 



which is equivalent to /(r+) = and corresponds to the situation when the two periods of the 
corresponding Weierstrass function become linearly dependent. This leads to a pair of solutions 

'"^^^ " -2he^^'t' ± (1 + (62 - 4ac)e2v^'^) ~ ' ^^^^ 



where 



6 = 4^3^, (68) 

c=2^-^-l. (69) 
' + 



3.2 5-D Schwarzschild null geodesies 

Here by 5-D it is meant 4 spatial dimensions. The relevant equation for u = 1 /r in this case is 

{u'f = 2Mu^-u^ + P, (70) 

where M is proportional to the five-dimensional mass. There is an interesting self-duality here, 
which is in fact a special case of Bohlin- Arnold duality, in that when we write the equation in 
terms of r we get 

{r'f =Pr* -r"^ + 2M . (71) 



10 



which is exactly the same with the constants interchanged. First consider the equation for u. 
Substitution of = ^Jjiu + ^) ^'^^^ take it into a form 

(y')' = 42/3 _ g,y - g, , (72) 

where 

9a = ^-8MP. (73) 

*^|(i-Mp), ,74) 

As usual, the r.h.s. of ([7^ has 3 real roots if 172 > and < (32/3)'^. The first condition is 
MP < 1/6 while the second is 8{MP)^ - {MPf < 0. Therefore we have 4 real roots if MP < 1/8. 
Also, note that the r.h.s. of equation (32) always has root —1/3 and expanding it into a power series 
around this point quickly shows that in fact 63 = —1/3. Finally, the point y = — 1/3 corresponds to 
r — 00. Hence in this case, with wi,W3 defined as before, we have 2 classes of solutions, depending 
on the initial conditions, scattering or trapped. The case of trapped paths is exactly the same as 
before, with the same addition formula for y and ?/((/)) = p{4> + wi). However, the scattering case 
is more interesting in 5D. This is because now the point in the C-plane corresponds to r = c» 
and so we can write the solution as y{(j)) — p{(p + UJ3) where the line </) = is in the direction of 
the ray incoming from 00 and G [0, 2u}i]. 

Things are even simpler when we solve the equation for r directly. The substitution = + 5) 
takes it into the equation (32) but now the difference is that point corresponds to r = 00 and so 
we can write the (scattering) solution simply as 



r(</.) = -^^p((/.) + i, (75) 



and the addition formula in this case is simply 



r(^i + 02) = -^Jf (piri^^W l,Pirih)y - ^ + J • (76) 



Finally, if MP > 1/8 then we have only 1 root of the r.h.s. of the Weierstrass equation and thus 
absorbing paths for which the solution is 

^ ^ ^p(0 + c^i) + iV (77) 



(r(0))2 2 V° 3 

where again the line (j) — is given by the direction of the ray incoming from 00. 

As before, we can get several special solutions by imposing — 27.g| which in this case gives 

MP = or MP =1/8. In the case MP = we get a special circular solution 

r{(j)) ^ V2M cos (j), (78) 

while in the case MP = 1/8 we get 

r(0) = 2VM(tanh(0/\/2))±i . (79) 

3.3 Duality and 7-D Schwarzschild null geodesies 

By Bohlin- Arnold duality [13], if we have a particle moving in Newtonian potential V oc r'^P~^ and 
following trajectory r{4>) = f{(j)) then there will be a particle with accordingly modified energy 

moving in a potential V oc r p following trajectory r(0) = /(0)^. In this case, ii V — —kr^P^^ 
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^ 2-2p ^ 

and particle has energy E then V — —Er p and E — k. 

Null geodesies in the (n + 1)D Schwarzschild geometry correspond to Newtonian motion in a 
r~" potential and so the duality applies to these geodesies as well. As we already saw, the 5-D 
corresponds to the case p = ~l and is self dual. A quick check reveals that the case p = —1/2 
gives a duality between null geodesies in 7-D and 4-D. 
Given the potential V = —kr~^, Newton's equation of motion is 

(r'f ^^r^-r^ + ^r*-". (80) 

The equation for null geodesies in 4-D is 

(r')2 = Pr-^ -r^ + 2Mr , (81) 

and in 7-D it is 

(r')2 = Pr'^ - + 2Mr-^ . (82) 

So, under the duality with p = —1/2 we have E ^ k and thus P -H> 2M. Therefore if we have 
4-D Black Hole with mass M and light with (E/L)'^ = P following the trajectory r{(j)) and 7-D 
Black Hole with mass P/2G and light with (E/L)'^ = 2M following the trajectory r — f{(j)) then 

2 



Making the substitution = y/P + 1/(3P) in the equation for r in the 7-D case takes it into the 

Weierstrass form with g2 — 4/(3AfP^) and 

53 = 8/27 - SMeP^. Therefore the orbits in 7-D satisfy 



r(0) = -^^p((/) + C) + i. (84) 



In this case, the formula for scattering paths looks especially simple, it is 



K0) = |yp(0) + ^- (85) 

By Bohlin-Arnold duality, the special solutions in 7-D corresponding to the special solutions in 
4-D given by P = 1/27 have 

P = -jL= (86) 

where M is proportional to the mass of 7-D black hole. The corresponding special solutions thus 
are 

r{^) = ^54M J i - \ . (87) 
y 3 1 ± cosh(0) 



3.4 Ellis Wormhole null geodesies 
3.4.1 Qualitative description 

The Ellis wormhole, is an ultra static solution of the Einstein equations coupled to a massless scalar 
field. While not necessarily physically very realistic, has been used in studies of gravitational lensing 
[5]. It has the metric 

ds^ = -dt^ + dr'^ + r(r - 2M){d9^ + sin^9d(t)'^) (88) 

Because qqq = —1, the physical spatial metric and the optical spatial metric coincide. Setting 
t = 0, 9 = ^ gives the optical metric on the equatorial plane. 
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If we set + y'^ = \/{r — MY ~ we may isometrically embed into E"^ with coordinates 
(a;, y, z) as the surface of revolution 



y^x^ ^ M cosh ^ , r = Af(l + sinh^). (89) 

Note that (j89p is a catenoid. This may be compared with the weU known Flamm paraboloid which 
gives an isometric embedding of the physical equatorial plane geometry of the Schwarzschild metric 



+ y2 = 2M + — , r=^x^+ 2/2 . (90) 

It is also possible to isometrically embed the Schwarzschild optical metric @ into Euclidean space 
but the formulae are more complicated: 



M,. A-r, 2M, 



y^^T^=^_=, 2^ / ^/_(4-9-)(l-— ) ^ (91) 



If we let u ~ then the equation of null geodesic is 



[u'f ^ (e - m^n^ + m + , (92) 

where ^ — Al"^ E"^ / . Note that this equation does not distinguish between r and 2M — r for 
r G [0, A/]. Before turning to the Weierstrass functions, we give a qualitative analysis of the null 
geodesies. Going back to the equation for r gives 

(^')' = -^r' - + (8^ - ly + 2M(1 - A^r + 2^/^(2^ - 1) ^ /(r) . (93) 
The roots of fir) have a very simple form, they are 

r = (l±i)M, (94) 

r^M{l± • (95) 

Extremal points of /(r), roots of f'{r) also have a simple form, they are 

r = M , (96) 



2^ 

From these result it follows that if 



r = M(l±W^-l). (97) 



^ S (0, 1/2) then /(r) has 1 positive real root M(l + -^/l/^ — 1) and 3 local extrema, all with 
value smaller than this root. 



• ^ g (1/2,1) then /(r) has 2 positive real roots M(l ± -^/l/C ~ 1) ^nd 1 global extremum 
(minimum) at r — Al. 

• ^ G (1, oo) then /(r) has no real roots and 1 global extremum (minimum) at r = Al. 
This shows that if 

• ^ G (0, 1/2) There are only scattering orbits with the distance of closest approach AI{1 + 
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• ^ G (1/2,1) There are both scattering and trapped orbits with the distance of closest 
approach M(l + -^/l/^ — 1) and the largest attained distance Af (1 — -^/l/^ — 1), respectively. 

• ^ € (1, oo) There are only absorbing orbits that is orbits incoming from oo that hit r = 0. 

There is an important point here. Suppose that we wanted to express r in terms of some Weierstrass 
function. The only way how to convert the full quartic into cubic is to substitute r = x + ro with rg 
being a root of /(r) — and then s — 1/x. If this approach is to be useful, we want rg G M, since 
otherwise, we would be looking for complex solution of the Weierstrass equation and the imaginary 
part C in p((j) + C) would not be half-period anymore but rather some analytically incalculable 
number and so this approach would not be useful at all. But /(r) has no real root in the case of 
absorbing paths and this foretells problems when treating this case. 

3.4.2 Weierstrass function approach 

First we make the substitution li^ = 1/a; in the equation (j92p which takes it into the form 

^ix' f = (e - l)x + (2^ - l)x' + -^x' . (98) 



Then the substitution 

takes it into Weierstrass form 
where 



M^y M2(1-20 



(99) 



{y'f - 42/3 _ g^y _ 53 , (100) 



ff2 = ^(l-e + e'), (101) 
53 = ^(2-3^-3^2 + 2^3). (102) 



Note that 52 > V$ and that 

(f)'-53^ = g(e-i)^e^>o, (103) 

unless ^ = 0, 1. Setting ^ = in eq. ([^^ shows that this case is not possible. The case ^ = 1 gives 
2 analytical solutions 

r.W-M(l±-^), (104) 

where comes from 00, r_ comes from r = and both are approaching r = M, but never reaching 
it. For other values of f the r.h.s. of equation (jlOOp has 3 real roots ei > 62 > 63 where 

ei =max(^^,^^^ , (105) 

e2 = mml— — - — 1, (106) 

63 = -^(1 + 0- (107) 
Now we will analyze the separate cases. Suppose that: 
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(0,1/2) . Then 

ei = ^ , (108) 

62 = , (109) 

63 = -^(1 + ?). (110) 

As a consistency check, one can verify that plugging y = ei into the expression r = r{y) 
indeed gives r = M(l + ^/^/^^^) as it should. Also, the point r = oo corresponds to the 
point y = oo and so the solution for the scattering orliits in this case is 



where the line ^ = is in the direction of the ray incoming from oo and </> e (0, 2u)i). Note 
that this solution always stays above r = 2M. 

The point r = corresponds to y = (5^ — l)/3 which is in this case in an unphysical region 
and so in accordance with Section 1.1 we only have scattering solutions in this case. 

i € (1/2,1) . Then 

61 = ^, (112) 
62=^, (113) 

63 = -^(1+0- (114) 

But now the scattering solutions penetrate into the region M < r < 2M and so I have to be 
careful here because r{y) is multivalued 



r = M(^l±-^^, + LJij . (115) 

This only becomes a problem once the orbit crosses r = 2M and so we didn't have to worry 
about it in the previous case ^ < 1/2. 

In this case (5^ — l)/3 > ei and r(5^ — l)/3) = or 2M. For orbits incoming from oo we 
clearly have to choose r(5^ — l)/3) = 2M because r = M is inaccessible. 
Also r(ei) = M{1 ± ^/l/£, — 1) and for the same reason we have to choose + for orbits 
incoming from oo. Hence as before 



where again the line = is in the direction of the ray incoming from oo and G (0, 2cji). 
What is left are orbits trapped in the region r < M(l — y/l/^ — 1). For these we need to 
choose minus signs in the above equations and so we get 



where now the additive constant in the argument of the Weierstrass function is necessary. 
This choice corresponds to setting the line </) = to be the axis of symmetry and 4> € (— /3, /3) 
where 

f°° dt 
= u>,- . (118) 

^(55-i)/3 \/ 4.1-^ - git - 93 
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^ g (l,oo) . Then 

ei = , (119) 

e2 = ^ , (120) 

63 = -^(1 + 0- (121) 

We know that in this case all orbits are incoming from oo and reach r = 0. Both y = 62 and 
y = 63 correspond to unphysical (complex) r and this time, y = ei corresponds to r = M 
without any ambiguity. Suppose we have an orbit starting at 00. (5^ — l)/3 > ei and so we 
need to choose plus sign in the relation r{y). 

Thus r{y = 00) = 00, then r{y = (5^ — l)/3) = 2M and finally we reach r{y = ei) = M. But 
if we continued the same Weierstrass function solution now r would begin to increase again, 
which we know is unphysical. Therefore we need to switch the branches and continue with 
minus sign in the relation r(y) so that we reach r{y = (5^ — l)/3) = 0. Now there is no way 
of continuing the solution and we need to start a new one, first using a minus sign and then 
a plus sign on its journey from r = to r = 00. Therefore an orbit going from 00 to r = 
travels a total angle loi + P and satisfies 



^ = l + -^yp(0) + ^ for <P€{0,oj,), (122) 



^ = l--^^p(</.) + ^ for </.G(a;i,/3), (123) 

where again line = is in the direction of the ray incoming from 00. 

General remarks 

(i) Note that the scattering solutions depend directly on p{(j)) and so the addition formula for 
Weierstrass functions can be applied directly. 

(ii) The same is true for the absorbing one, however we need to be careful to stay in the region 
(j) € (0,u;i) or ^ € (wi,/3) when applying it. 

3.4.3 Angle of deflection in the scattering case 

The equation for u can be factorized as 

{u'r = (1 + M^u') + (^ - 1)^^') • (124) 



Now, the distance of closest approach is ro = M + My/l/^ — 1, which corresponds to 

(125) 



(126) 




My 1 

Let / be half of the angle (f> travelled by the light. 

du 



Jo 



Vl + M^uVC/^^ + (C - ' 
Making the substitution u = uot, we have 

'■1 dt /"^ dt 




(127) 
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Write for the final integrand above. It is straightforward to differentiate / n times w.r.t. f 

and the resuh is 

-J- = (\-e) — \ — — I , 128 



Therefore we can expand / as 

/(i,0 = f:(l-t')'^^^^^r^2-"r. (129) 

n=Q 

Scattering orbits exist for ^ e (0,1) and for this range of values of f the sum converges uniformly 
(for example by straightforward application of the Weierstrass M-test) and therefore we can write 

I^y( illl_m2-^C f\l - t'f-^] . (130) 

The integral in this sum can be computed by hand, one way is as follows. The volume Vn of an 
n-dimensional ball 

Vn = —, r (131) 

r(i + f) 

Therefore 

^^"^^ ~2F2„_i~2r(n + l) ^"-1/2 -2n!^(v 2;' ^ ^ 

Now using the identity 

r = (2n- l)!!2-"Vi (133) 

gives 



ri=0 ^ ' 

This can be further simplified using the identity (2n — l)!!n! — (2n)!2~" to give 



^ = ; E . (135) 



ri=0 ^ ^ 

Now the angle of deflection S(j> is given hy S(f> = tt — 21 and so 

/ 2n^ 



^^ = 7r-^^ I ^ j 2-4"^- (136) 

The first few terms of this expansion are 

Scb = --f - — f2 _ _ 1225^,4 _ 396971 5 _ 53361^ ^ _ .^37^ 

^ 4^ 64^ 256^ 16384^ 65536 1048576^ ^ ' 

with C = (M/6)2. 

We have also tried expanding the defiection angle in terms of /i M/tq following [B]. Substituting 

e = 7^ 72 (138) 

i -1 
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into the integral (I127p and expanding in the powers of ^, using Mathematica, the first few terms 
are 

(50 1 2 1 3 41 4 9 . 25 6 37 7 11959 = 1591 g 13311 ^ 29477 ^ 

— = U M M M H M H M H A* H M -■ • ■ ■ 

TT 4^ 2^ 64^ 16^ 256^ 128^ 16384^ 2048^^ 65536^ 32768'^ 

(139) 

This expansion is not very useful , since the coefficients don't seem to be decreasing very fast, the 
coeff. of /i^^ is almost 1. 

Note that this expansion is completely different from the one given in [B]. 

4 Conclusion 

In this paper, we have used Weierstrass elliptic functions to give a full description and and classifi- 
cation of null geodesies in Schwarzschild spacetime. We then used this description to derive some 
analytical formulae connecting three points on these geodesies and found second order expansion 
of the deflection angle in the scattering case. Finally, we derived some properties of light triangles 
in this spacetime and used the Gauss-Bonnet theorem to derive a quantity which gives the same 
answer when integrated along a scattering geodesic, independently of the geodesic in question. 
We then showed that the Weierstrass elliptic function formalism can also be used to describe other 
more exotic spacetimes such as Reissner-Nordstr0m null geodesies and Schwarzschild null geodesies 
in spacetimes with spatial dimensions 4 and 6. In all these cases, the elliptic function approach 
allows one to find the special case analytical solutions with ease (simply by looking at the values 
of parameters for which the elliptic function in question collapses into a periodic one). 
Finally we applied the formalism to describe the null geodesies of the Ellis wormhole and found 
an expansion for the angle of deflection in this case. 

After the appearance of the first verson of this paper on the archive, Betti Hartmann pointed 
out to us that our results may easily extended to the case of a Schwarzschild black hole pierced 
by an infinitely cosmic string studied in |15) . One need only replace the variable by 50 , where 
< (5 < 1 is the deficit parameter. Similar remarks to the other metrics studied in this paper. 
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